NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



IVIASA 

Technical Memorandum 81985 

(tiASA-TA-81985) JUA SI-Si.rtUl.TA NEOUS 880-31285 

0B3EHV AX IONS OF THE BL LdC OBJECT M 501 It) 
j BAY, UV, VISIBLE, IB A80 BAOIO FHEQU EEC IBS 
' KSA) 32 p HC A03/Mt A 01 CSCL J3A Un-.Ias 

G3/89 31330 

Quasi-Simultaneous 
Observations of the BL Lac 
Object MK 501 in X-Ray, UV, 
Visible, IR and Radio Frequencies 

Y. Kondo, D. M. Worrall, R. F. Mushotzky, 
R. L. Hackney, K. H. Hackney, J. B. Oke, 
H. Yee, G. Neugebauer, K. Matthews, 

P. A. Feldman and R. L. Brown 


JULY 1980 


National Aeronautics and 
Space Administration 

Goddard Space Right Center 

Greenbelt, Maryland 20771 



QUASI -SIMULTANEOUS OBSERVATIONS OF THE BL LAC OBJECT MK 501 
IN X-RAY, UV, VISIBLE, IR AND RADIO FREQUENCIES 


Y. Kondo 

Laboratory for Astronomy and Solar Physics 
and 

D.M. Worrall^ and R.F. Mushotzky 

Laboratory for High Energy Astrophysics 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 

R.L. Hackney and K.R.H. Hackney 

Department of Physics & Astronomy 
Western Kentucky University 

J.B. Oke and H. Lee 
G. Neugebauer and K. Matthews 

2 

Palomar Observatory 

California Institute of Technology 

P.A. Feldman 

Herzberg Institute of Astrophysics 
NRC of Canada, Ottawa 

Robert L. Brown 

3 

National Radio Astronomy Observatory 
Green Bank 


^Also University of Maryland. Now at University of California, San Diego 
2 

Observations were also made at the Mt. Wilson Observatory as part of a 
collaboration between California Instiwte of Technology and Carnegie 
Institution of Washington. 

3 

Operated by Associated Universities, Inc., under contract with the 
National Science Foundation 


ABSTRACT 


Quasi -simultaneous observations of the BL Lac object MK 501 were performed 
for the first time at ^-ray, ultraviolet, visible, Infrared and radio frequencies. 

As the BL Lac objects are known to vary In their flux, such a "quasi -instantaneous" 
spectral erergy profile Is necessary In order to describe properly the energy 
generation mechanism. Th'- observed spectral slope from the X-ray to UV regions 
Is positive and continuous, but that from the mid-UV to visible light region 
becomes gradually flat and possibly turns down toward lower frequencies; the 
optical-radio emission can not be accounted for by a single power law. Several 
theoretical models have been considered for the emission mechanism. In some 
cases quantitative comparison with the data Is not practical. However, most 
of the models are, at least, not inconsistent with the observations. A quantitative 
comparison has been performed with the synchroton-sel f-Compton model; the total 
spectrum is found consistent with this model. The spectrum from visible light 
to X-ray is consistent with synchrotron radiation or with Inverse-Compton scat- 
tering by a hot thermal cloud of electrons. The continuity of the spectral 
slope from X-ray to UV Implied by the current data suggests that the previous 
estimates of the total luminosity of this BL Lac object has been underestimated 
by a factor of about three or four. 


I. INTRODUCTION 


BL Lacertae objects share a number of properties In common with quasars 
and, as a consequence, are generally considered to be closely related to quasars. 

Features of BL Lac objects by which they are often distinguished from "average 

quasars" are as follows. Many 8L Lac objects have been located in the cores 

of elliptical galaxies; the spectrum of the core region exhibits little or no 

detectable line emission indicating a virtual absence of gas; the optical and 

radio emissions are highly linearly polarized; the fluxes vary significantly 

at all wavelengths. 

Relatively larger and more rapid changes have been observed in visible 
light and in the near-infrared regions than at the radio wavelengths. In visible 
light, variations have been observed on a time scale on the order of 10 days 
(Stein, O'Dell and Strittmatter 1976). Based on archival plates, there are 
apparently quiescent and active phases in some sources, which may brighten 
by as much as some 5 magnitudes during the active phase.. Variability in the 
ultraviolet and X-ray regions is currently being investigated using the satellite 
observatories. In the case of MK 421, the 2-10 keV X-ray flux was observed to 
change by an order of magnitude on the time scale of days (Ricketts, Cooke and 
Pounds 1976). Modulation collimator experiments have recently verified and 
established unambiguously the following five BL Lac objects as X-ray sources: 

MK 501, MK 421, PKS 0548-322, PKS 2155-304 and 2A 1219+305 (Hearn, Marshall and 
Jernigan 1979; Schwartz et al . 1978, 1979). 

It is generally suspected that the energy generation mechanism in quasars 
is the same as that in BL Lac objects. Since the spectrum of a BL Lac object 
is relatively uncontaminated by the spectrum of gas, the physical mechanism 

for the non-thermal radiation in both types of source may be investigated 
effectively by observing BL Lac objects. These observations must be made on a 
timescale short compared to the characteristic timescale for variation of the 
source. Observations of both the flux and polarization of 1? BL Lac objects by 
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Angel et al . (1978) suggest that the shortest characteristic timescale for the 
erratic variability seen In some objects is £ 1 day. Similarly, In the X-ray, 
the shortest observed timescales for a BL Lac object, PKS 2155-304 (Snyder et al. 
1980) is approximately one day. However such rapid variability has not been 
detected In MK 501. For this object no variability has been reported at radio 
frequencies and the variability at visual wavelengths reported by Barbierl and 
Romano (1977) occurred on a timescale of years with a range of .9 mag. 

Our present observations Indicate that the X-ray variability also occurs 
on timescales longer then one week since no significant intensity change was 
seen during each of our 6 day observation periods while significant spectral 
changes were seen over a 6 month timescale. We therefore consider the radio, 
Infrared, optical, ultraviolet and X-ray observations reported In this paper to 
represent a quasi -simultaneous "snapshot" of the total spectrum of MK 501. 

II. OBSERVATIONS 

The data discussed in this work were obtained in August and September 
1978 with the HEAO-1 and IUE satellites, the 5-meter and 2.5-meter telescopes 
at the Hale Observatories, the NRAO 92-meter radio telescaope at Green Bank, 
and the 46-meter radio telescope at the Algonquin Observatory. 

The X-ray measurements were made with the A2 experiment on the HEAO-1 
spacecraft. The detectors are comprised of mechanically collimated gas-filled 
multi-wire proportional counters (see Rothschild et al . 1979 for the experiment 
description). Measurements were made with the argon and xenon-filled counters. 
For a spectral determination in the 2-50 keV range, MK 501 was observed in a 
pointing maneuver of the satellite between 1978 September 8 19:00 and September 
9 01:00 UT corresponding to JD 2443760.29 and JD 2443760.54. Net source counts 
were obtained by subtracting background counts accumulated from a region close 
to the source which the detectors w°re scanning between 1978 August 21 and 
August 31 UT. During this scanning period, MK 501 was in view for about 1 minute 
every 30 minutes, enabling its intensity to be monitored throughout the interval. 

x , K - — * 
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The IUE Instrumentation has been described by Boggess et al. (1978). 

For this program, spectral observations of MK 501 were made In the low-dispersion 

o 

mode (effective resolution,^ 10 A) using the optimally configured SWP (short- 
wave length prime) and LWR (long-wavelength redundant) cameras. The total 
frequency (v) range observed was 14.97 < log v «. 15.36 MK 501 was observed 
using the 10 x 20 arc sec aperture to insure acquisition by blind offset and 
to permit application of the flux calibration which is appropriate for this focal- 
plane aperture. The cameras were prepared using both the XPREP sequence (to 
erase all traces of previous over-exposures) and the SPRFP sequence (to minimize 

the unexposed background level). The short-wavelength spectral image SWP 2394, 

o 

covering the range between 127b and 1910 A, was obtained 1978 August 25 03:22 - 

o 

09:42 UT; the region shortward of 1275 A was contaminated by geocoronal and 

interplanetary Lyman -a emission. The long-wavelength spectral image LWR 2177, 

o 

ranging from 1900 to 3140 A, was obtained 1978 August 26 03:12 - 08:12 UT 

Both IUE exposures were made during the lowest-radiation part of the IUE orbit. 

Under these conditions, the radiation-contributed background over the relatively 

long exposure times amounted to approximately one-quarter the total response 

at the wavelengths of peak sensitivity of the cameras. The effective 

o o 

signa 1 *:o-noise ratio for single data elements (2 A bandwidth SWP, 4 A bandwidth 

LWR) indicated by the standard deviation of the measurements of net signal after 

subtraction of the smoothed background is typically about 12:1. 

The 5-meter Hale telescope observations were made using the multichannel 

spectrometer for visible light with an aperture size of 10 arc seconds. The 

o 

frequency range covered was 14.42 < log v < 14.96, with 40 A bandpasses for a < 
oo o 

5700 A and 80 A for A > 5700 A. Similar observations were obtained on the following 
dates: 1978 May 26 9:25; 1978 Aug,, 29 4:15; 1978 Aug., 30 3:39; and 1978 Aug,. 

31 3:34 UT. 
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The Infrared observation at Mt. Wilson was made with the 2.5-meter 
telescope at 2.2 urn (ax * 0.4 urn) on 1978 August 31 03:08 UT. The aperture was 
15 arc seconds in diameter. 

Observations with the 46-meter telescope of the Algonquin Radio Observatory* 

* 

The Algonquin Radio Observatory is operated by the National Research Council of 
Canada, Ottawa, as a national radio astronomy facility. 

were made at 10275 MHz (log v » 10.02) (1978 August 15-24 UT) and 10650 Wz 
(log v * 10.03) (1978 August 30 - September 3 UT), using a beam-switched receiver 
with system temperature of about 150 K and bandwidth of approximately 100 MHz. 

The 92-meter radio telescope observations at Green Bank were obtained on 
most days in the period 1978 August 15-31 at source transit (between 12:00 and 
13:00 UT). The frequency used was 4750 MHz (log v * 9.68); the system temper- 
ature was 70 K and the instantaneous continuum bandwidth was 500 MHz. 

III. REDUCTION AND DESCRIPTION OF THE DATA 
A. X-Ray Spectrum (HEAO-1) 

The X-ray spectrum was derived from the HEA0-A2 pointing observation 
between 1978 September 8 19:00 and September 9 01:00 UT. The spectrum is soft 
enough for the best statistical accuracy to be achieved with the argon detector. 
The upper portion of Figure 1 shows the incident counts along with the best model 
fit. Xenon detector data give consistent spectral fits. Figure 1 also shows 
the inferred incident spectrum after folding through the detector response. The 
best value for the index of the differential photon spectrum, r, is ^ 2.52 and 

21 

the best fit low energy absorption implies a column density, N^, of ^ 4.6 x 10 

_2 

atoms cm . In this paper we use a for the slope of a power law fit to the 

continuum such that the flux at a given frequency f is defined as f = Av" a 
-2 -1 

(Wm Hz ). The photon index r, is related to the energy index a by r = a + 1 . 

_ r -2 -1 -1 

The X-ray spectra are fitted to a form f = Av exp(-aN^) ph cm sec keV 
where is the equivalent column density of hydrogen atoms assuming that the 
material is neutral and 0 is the appropriate cross section (Fireman 1974). The 


90% confidence error contour for r and is also shown In the figure. No 
additional structure to the spectral form significantly Improves x • A thermal 

O 

bremsstrahlung fit Is found to give a worse x • although one which Is acceptable 
at 70% confidence limit. The best temperature In this case Is 4.5 keV. 

In the HEAO-1 nominal scanning mode, the spin axis points towards the 
sun, thus moving a degree each day, and the A2 detectors take ^ 30 minutes 
to trace each great circle on the sky. MK 501 Is at fairly high ecliptic latitude 
and so was In view for about one minute during each scan for 10 days (1978 August 
21-31). The courage Is very Incomplete because of pointing operations at other 
regions of the sky for several periods of 8-12 hours during this time. There 
are too few counts to enable spectral parameters to be determined with sufficient 
precision to discern any spectral change between this period and that In September. 

For the pointing observation, the 2-6 keV and 2-20 keV energy fluxes are 
(2.9^°) x 10” 11 erg cnf 2 s _1 and (4„5* 2,, g) x 10" 11 erg cm" 2 s‘* 1 , respectively. 
Experimental systematic effects and the unavoidable presence of Intrinsic sky 
fluctuations dominate the quoted errors. The average intensity from scanning 
data can be determined with comparable precision and Is found to be in agreement. 
Furthermore, one-quarter day averages for the August scanning observations are 
consistent with a constant flux to within 60% confidence. 

If we compare our observations with previous X-ray observations of MK 
501, we find no evidence so far for changes in the 2-6 keV energy flux. The 
revised UHURU value of (2.75 + .52) x 10"^ erg cm" 2 s"^ (Forman et al . 1978), the August 
1977 HEA0-A3 value of (3.6 +1.2) x 10"^ erg cnf 2 s"^ (Schwartz et al. 1978) 
and the March 1975 Ariel 5 value of (3.6 +1.8) x 10“^ erg cm" 2 s - ^ (Snijders 
et al, 1979) are comparable to the present value. Simultaneously with the A3 
measurement, the A2 detectors observed the source in August 1977 (Mushotzky et 
al. 1978; Mason et al. 1979). The 2-6 keV flux agreed with the A3 value and 
hence with the present one. However, the spectrum was of a different form. 

A single value of r, as defined earlier, did not give an acceptable fit to the 
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data. A photon Index r - 25 from 0.15 to 3 keV, with a hardening to about 1,2 
above 3 keV, was required. Mason et al, (1979) suggest that this two component 
nature may also be present In MK 421 , They Interpret the large spectral and 
Intensity change in MK 421 which occurred over 6 months (Mushotzky et al. 

1979) as due to variability of the hard component by about a factor of 5, 
combined with a relatively constant low energy component which extends from 2 
keV to lower energies The present MK 501 observation suggests that this object 
also may always be characterized by two components, where both are seen to change 
in intensity The suggestion is that between August 1977 and August 1978 the 
hard power law component with r ~ 1.2 decreased by about a factor of two in 
intensity while the soft component with r 2 5 increased by a similar factor. 

The 1.2 - 12 keV Ariel 5 data were shown by Snijders et al. (1979) to be incon- 
sistent at 90% confidence with single power laws of slope less than 1 2 or greater 
than 2 3, and gave r = 1.8 + 0 5 as the best fit However, these authors have 
not considered a two -power law fit Inspection of their data would not appear 
to rule out such a possibility, in which variability of the two components would 
be suggested. While it is premature to draw strong conclusions, we tentatively 
suggest that BL Lacertae objects as a class are characterized by such multiple 
component X-ray spectra, where the overlap is seen at a few keV In contrast, 
when the 2-50 keV flux from the quasar 3C 273 was examined for spectral changes 
accompanying the intensity changes discerned over a few months, none were found 
at a confidence level comparable to those reported here for spectral changes 
in MK 501 (Worrall et al. 1979) 

8. Ultraviolet Spectrum (IUE) 

The IUE data for MK 501 were calibrated to absolute fluxes using the 
calibration established by Bohlin et al (1979), Additionally, special obser- 
vations of calibration stars were required for deriving corrections specific 
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to the MK 501 exposure conditions to remove effects Introduced by on erroneous 
photometric Intensity transfer function used for processing SWP Images prior 
to 1979 July. The calibrations and the Image processing procedures are continually 
subject to evaluation and possible revision as the observing experience with 
the IUE Instrumentation progresses. Uncertainties In the calibrations are cur- 
rently estimated to be no greater than about + 15%. 

Calibrated data from IUE exposures SWP 2394 and LWR 2177 are listed In 
Table 1 and plotted In figure 2. Individual data samples have been binned into 
larger wavelength Intervals to increase the effective signal -to-noise ratio 

for determining local continuum levels Included measurements cover ranges of 
o o 

3140 - 1900 A (LWR camera) and 1910 - 1275 A (SWP camera). No discrete spectral 

features of the source were evident In the fully resolved data set, for which 

o 

the detection limit is estimated to be ^ 4 A in equivalent width. 

A heavily exposed IUE observation of the BL Lac object MK 421 which, like 

MK 501, has a detectable galactic envelope In visible wavelengths, revealed no 

evidence of flux from the surrounding galaxy in the wavelength range of the 

o 

LWR camera (3200 - 1900 A). Thus, we feel that the IUE spectra of MK 501 are 
essentially due only to the compact nuclear BL Lac source with little contribution 
from the associated galaxy. 

In Figure 2 are plotted the IUE data both before and after the application 
of a correction for effects of interstellar reddening. Since no direct deter- 
mination of reddening is available for MK 501, the mean color excess E(B-V) = 

0.083 mag appropriate to its galactic latitude was adopted as the basis for 

o 

applying the ultraviolet extinction curve of Code et al. (1976). The 2170 A 
extinction feature is not significantly detectable in our observation with a 
substantially better signal-to-noise ratio than that of the short-exposure data 
used by Snijders et al. (1979). The weakness of this feature suggests that E(B- 
V) is less than 0.1 mag in the direction of MK 501. Snijders et al. adopted 
a larger value for the color excess in order to minimize an apparent depression 
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near the center of tholr IUE spectrum. For our calibrated and dereddened 

spectra, we obtain the following power-law Indices (In the sense f « v’ a ): 

o 

LWR a * 0 31 ♦ 0 23 (3140 - 1900 A) 

SWP a « 1.23 + 0.35 (1910 - 1275 A) 

o 

LWR ♦ SWP a » 1 02 + 0 14 (3140 - 1275 A) 

The quoted errors are 95% confidence limits for the power-law Indices estimated 
on the basis of the internal error of the data sets; they do not include the 
uncertainties in either the reddening correction or the calibrations. If the 
appropriate E(B-V) is actually less than the adopted value of 0 083 mag, then 
the ultraviolet continuum of the source is steeper than indicated by the results 
here reported The degree of possible reddening of the BL Lac source due to 
its surrounding galaxy is also unknown. The significant difference between the 
indices of power laws fitted to the separate LWR and SWP data sets is interpreted 
as being due to a curvature of the spectral energy distribution, with a turnover 
toward longer wavelengths. These results differ significantly from the 
a = 0.50^Q - ]g result of Snijders et al . (1979) for the total UV spectrum. 

Ci Visible Light Spectrum (5-Meter Telescope) 

All the observations give the same magnitude (V ~ 13,75) and energy dis- 
tribution, A typical spectral scan is shown in Figure 3, The spectrum clearly 

contains a galaxy, as one can see the H and K Ca II break, the depressions at 
o o 

4150 A, 4300 A, the Mg-b band, and effects of bands in the red, A match up with 
a standard elliptical galaxy yields the redshift z = 0 034 + 0,002, The energy 
distribution is flatter than a normal elliptical and the H and K lin^s are ab- 
normally weak, as is typical of galaxies with a nonthermal blue continuum, 


- 11 - 


A model fitting procedure was performed by adding a normal elliptical 
<r*l "<y energy distribution (Yee and Oke 1978) to a power law of Index a. There 
are no reasonable fits when a Is +1.0 or +0.5. The best fit occurs with a • 0; 
but a ■ -0.5 Is also possible. A very small reddening correction corresponding 
to Eg_ v ■ 0.05 Is applied. 

Fits are all given In terms of visual magnitude and for a ■ 0. For the 

10 arc second aperture the non-thermal part has m y * 15.11, which corresponds 

to 0.33 y ID" 28 W m’ 2 Hz'*. The galaxy part has m y * 13.94 or 0.96 x 10" 28 
-2 -1 -1 -1 

W m Hz . For H Q ■ 55 km s Mpo and q Q * 0, the absolute visual magnitude 
of the galaxy Is -22.3 within a 10 arcsecond diameter aperture. The galaxy 
appears to be highly concentrated towards the center since the magnitude dif- 
ference for a 7" and 10" aperture is only 0.18 mag. Thus the total absolute 
visual magnitude of the galaxy Is perhaps -22.5 to -22.8; the galaxy Is therefore 
of only modest luminosity. 

Previous measurements of the continuum by Maza, Martin and Angel (1978) 
indicate that the polarized continuum radiation has a slope a = 0.83 + 0.10, 
considerably steeper than the fit to the continuum derived here. This fit is similar 
to the overall fit to our UV data. This suggests that the non-thermal component has 
varied between the epoch of Maza et al. and the present observations. Compari- 
sons of fluxes in the V, B ar.d U bands (log v ^ 14.73, 14.83, and 14.92 respect- 
ively) show a possible change in the B flux from Maza et al.'s value of log ~ 
-28.02 W m’ 2 Hz"^ to our value of log F^ = -28.19 (see Figure 5 of Maza et al . and 
Figure 3 of this paper), while the V and U fluxes were the same. We also note 
that Maza et al.'s values of U- B-V, and U-B differ from those of Tapia, Craine, 
and Johnson (1976) thus giving furtt’ " evidence for optical intensity and 
spectral variability. 
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D. Infrared Data (2.5 Meter Telescope) 

The procedures for the data reduction are as described by Neugebauer et 
al. (1978). 

The measurement at 2.2 urn at 15" diameter aperture yielded a flux value 
of 35 + 3 m Jy or ■ 10.6 ♦ 0.1 where the error Indicated Is la. If the 
galaxy part of the visual light Is extrapolated to a 15" aperture using the 
same curve of growth that exlits between 7" and 10", then m / ■ 13.75. This 
produces a color In close agreement with the average value of V-K ■ 3.3 for 
elliptical galaxies found by Frogel et al . (1978). If allowance Is made for 
the uncertainties In the photometry and the extrapolation to m y and, also. If 
V-K ■ 3.0 Is taken as a lower limit for elliptical galaxies, then roughly one 
half of the 2.2 pm flux In a 15" aperture could be non-thermul . This produces 
a lower limit of a < 1.2 between 0.5 and 2.2 um. This limit on a Is c r, sistent 
with the value of a derived In this paper from the visual data and a * .82 
between 0.44 and 10.6 pm found by Ulrich et al . (1975). 

E. Radio Data (46-Meter Telescope) 

Throughout the period of the observations pointing corrections were 
checked, gain stability was monitored, and focus offsets were maintained by 
measurements of strong point calibration sources.. Ar.y residual pointing errors 
should not have exceeded about 20", which for the 2 '40" beamwidth, corresponds 
to ^ 5% of the measured flux density 

The observational n .'suits are plotted in Figure 5 Over the time inter- 
val 1978 August 15 - September 3 UT, there is no evidence for radio vari'vbil i ty 
of MK 501, in agreement with the 92-meter observations at 4750 MHz. The mean 
flux densities are 1 14 + 0.02 Jy at 10,275 MHz, 1.14 + 0.01 Jy at 10,650 MHz, 
and 1.14+ 0.01 Jy overall at 2 85 cm, where the errors quoted are at 2a These 
results are identical, within experimental error, to earlier measurements at 
28 cm (see, e.g,, McCutcheon and Gregory 1978) 
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F, Radio Data (92-Meter Telescope) 

The results of the radio data at 4750 MHz are plotted in Figure 4. The 
error bars are + 2a. Over the time interval of these observations the source 
was quiescent,, The mean flux density is 1 ,327 Jy, a value in keeping with 
the value of 1.295 Jy interpolated at this frequency from the observations 
at 2.6 and 8.1 GHz made by Sramek and Toumassian (1976). 

IV. DISCUSSION OF OBSERVATIONAL RESULTS 
Constraints from the Overall Spectrum 

We have plotted all of the quasi -simul taneous data for MK ^01 in Figure 
6 As was discussed earlier, it is probably safe to assume that the contribution 
of photons from the galaxy is negligible to the X-ray and ultraviolet data. 

In visible light there is a substantial contribution from the galaxy. We have 
therefore plotted both the observed flux (galaxy and the BL Lac source combined) 
and the deduced BL Lac source flux; the forme” is indicated by a line and the 
latter by a broad bar. No correction for the galaxy light can be made for the 
IR data at the moment. The radio flux may be considered to be basically all 
from the BL Lac source. 

The most significant feature in the observed spectral energy distribution 
is the turc.-'ver from the UV into the visible light component. We should also 
note that although there exists an apparent continuation of the slope from the 
UV to the X-ray region, the slope for the intervening region (a 0.9) fits the 
UV data but not that within the observed X-ray band. We infer a steepening of 
the slope at X-ray wavelengths, which is possibly related to the electron energy 
loss mechanism in the source. 

Various emission mechanisms have been proposed for the BL Lac sources, 
including synchrotron self-Compton emission (Jones, O'Dell and Stein 1974b), 
Compton emission (scattering) from thermal electrons (Katz 1976), bremsstrahlung 
emission from accretion disks (Eardley et al. 1978), and emission from turbulent 
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plasmas (Colgate et al. 1975) or hot gas clouds (Fabian et al. 1976; ^chnopper 
et al. 1977). We shall concentrate on the first two mechanisms because they 
have been presented In forms amenable to quantitative assessment. 

A. Synchrotron Self-Compton Emission 

Two features of the total flux density spectrum (Fig. 6) are Immediately 
suggestive of a synchrotron self-Compton model for the radio and visible-light, 
ultraviolet. X-ray emission from MK 501: these are (1) the flattening or turn- 

over In the radio spectrum at v * 10 GHz [see Schwartz et al. (1978) for a more 

complete representation of the radio spectrum J and (2) the corresponding flat- 

1 5 

tening in the visible light at v s 10 Hz. This Is precisely the shape of the 
total synchrotron self-Compton spectrum expected from a source in which the radio 
emission is synchrotron radiation while the visible light, ultraviolet and 

X-ray emission (all frequencies ^ 10 14 - 10^ 8 Hz) is first-order Compton scattering 
of the relativistic electrons with the radio photons (Jones et al . 1974a, 1974b; 
Burbidge et al. 1974). Since definite relationships exist between the synchrotron 
and Compton components we can readily Investigate the consistency of this model 
in the particular case of MK 501. 

The first-order Compton flux density can be expressed in terms of the 
observable radio properties of the assumed spherical, isotropic and homogenous 
source (Marscher et al. 1979) as 


S c (v) 


10" 4 K s (a) In 


v _2 - 2(a+2) -2(a+2) u -(3a+5) 

v m a m 0 v m 
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L 


2.4 x 10 
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J 


10- 26 W cm~ 2 s"^Hz‘ 


where is a function tabulated by Marscher et al. (1979) and is of order 1. 

Here, is the upper frequency limit of the synchrotron radio spectrum, 

v m (GHz) is the frequency where the radio flux density is maximum, S m (Jy) is 
the radio flux density at u , e (mill iarcsec) is the angular diametsr of the 
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source and v (Hz) Is the frequency of the first-order Compton scattered photons. 

A 

For a consideration of the effects of relativistic motion and non -homogen Ity, 
see Gould (1979) and Jones, O'Dell and Stein (1974a). For the rest of this 
discussion we shall assume such effects are small. The spectral Index used In 
this expression, a Is the spectral Index of the optically thin radio spectrum 
which is difficult to establish because the source appears at least partially 
optically thick at all frequencies less than 10 GHz (l.e., the radio spectrum 
Is flat at these frequencies). However, In the synchrotron sel f-Compt.on model, 
the optically thin synchrotron spectral Index is identical to the Compton spectral 
Index and this latter quantity we can establish from our (Fig. 6) UV and X-ray 
flux densities as a s 0.9. This value Is different from the local index in either 
the UV or the X-ray band and is appropriate for an UV-X-ray Index, But, our 
result Is insensitive to the exact values of the spectral index «. Using this 
value of the spectral Index together with the angular diameter of the radio core 
of MK 501 as measured by VLBI techniques at 6 cm, e = 1.1 mas (Weiler and Johnston 
1980), R ^3.4 l.y., and the extrapolated maximum flux density S m s. 2 Jy 
at a turnover frequency of the core, v m s 2 GHz, we obtain an X-ray flux of 

S r (lu 18 Hz) a 2.2 x 10 32 W m“ 2 Hz -1 , 
which is in good agreement with the observations. The magnetic field strength 
Is also determined in this model, 

B = 3.4 x 10‘ 5 e 4 v m 5 S m " 2 s 4.0 x 10" 4 G. 

The flattening of the first-order Compton emission that appears at visible 
light wavelengths can be accounted for as a low energy break or turnover in 
the Compton spectrum. Briefly, if the power law distribution of relativistic 
electrons has a low energy cutoff at an energy , then the self-Compton radiation 
spectrum will reflect this cutoff at a frequency 


= y 


4 E, 
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For MK 501 if v ' s 10* 5 Hz, then Ej s 310 MeV. These electrons will radiate 
synchrotron emission at 200 MHz and are therefore not — esslble to direct obser- 
vation In the raiilu spectrum (see Marscher f 1 977 ) for a description of the detailed 
spectrum for this case). The first < rue 1 ' 'ompton radiation will extend from 
^ 10 15 Hz to an upper frequency limit v £ , which *s determint-d by the high energy 
cutoff of the distribution of relativistic electrons. Let this high energy limit 
be E 2 » then, 

E ' 2 

2 4 2 , 

v c 2 T "T v 2° 
me I 

Choosing Vg conservatively as the highest frequency at which radio emission 
has been detected from MK 501, Vg = 100 GHz (Joyce and Simon 1976) 
we find that Eg = 4 GeV (from the synchrotron critical frequency in a 

field of 4 x 10" 4 G) and hence that v £ 2 s 8 x 10^® Hz. At frequencies near v c 2 

we expect the spectral index to begin to steepen; the X-ray data are consistent 
with this predicted steepening. 

Finally, the total self-Compton luminosity 10 15 - 10 19 Hz is L s 5.5 

45 -1 C 

x 10 ergs s assuming that the distance to MK 501 is ^ 200 Mpc as given by 

the redshiftfe = 0.0337 (Ulrich et al. 1975) and z = 0.036+0.002 (present work)), with 

H Q = 55 km s” 1 Mpc"^. The energy density in relativistic electrons is 

u p 3 1.6 x 10" 2 e" 9 v m - 7 S m 4 = 8.2 x 10"^ erg cnT^, 

which means that the total energy in relativistic electrons is U g « 1.2 x 10 94 

ergs and the resultant lifetime of the radiating particles t = U g /L e = 6.8 

years. This timescale compares favorably with the timescale for X-ray 

variability (Mushotzky et al. 1979). 
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The lifetime for spectral variability in the X-ray can be estimated 

from the effects of synchrotron and Compton losses on the electron spectrum. 

The lifetime for synchrotron emitting electrons whose Compton scattered photons 

are observed at 5 keV is ^ 4 x 10 years, while the Compton loss lifetime is 

^ 4 years,, The high frequency break in the photon spectrum caused by Compton 

IT-? 

losses appears at \>g = 3,6 x 10 t Hz, where t is the time after injection 
in years. So, the spectrum is expected to break in the far IR region. As the 
preceding equation shows, the lifetime is a very sensitive function of o and 
v and the fact that the calculated lifetimes approximately agree with the obser- 
vations is somewhat fortuitous. 

We note that our solutions differ significantly from the similar calculations 
of Schwartz et al„(1978). This is primarily due to the recent angular size data of 
Weiler and Johnson (1980) and reflects the problems inherent in deconvolving 
radio spectra composed of multiple components., 

B„ Inverse-Compton Scattering Off Thermal Electrons 

The apparent continuation of the power law between the visible light and 
X-ray regions suggests that the mechanism discussed by Katz (1976) may be appropriate 
for the source. In his mode 1 sion is due to Compton scattering. The origin 

of the initial supply of soft photons is not explained. The electrons originate 
in a hot gas cloud and have non-relativistic velocities. Radiation is restricted 
to below a few hundred keV since the Comptonized photons must have energy less 
than that of the non-relativistic electron gas. A fit to the observed spectral 
slope, allowing variability by as much as a 25% change in flux in a period over 

O „ '3 

the order of one year, gives an electron density 10 cm and a source size 
0.3 light years. The mass necessary to bind the gas cloud at the limit of maximum 
source size is 3-10 6 M .. 

C. Synchrotron Radiation 

The appearance of a "break" by approximately one unit in the visible=l ight- 
UV spectrum and the extension of that "broken" power law into the X-ray band 
suggests synchrotron radiation as the source of the visible*! ight-UV-X-ray 
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continuum. In this model the radio flux would be from a different, presumably 
larger component. If we regard the "break" as being due to synchrotron energy 
losses and use •v 1.5 x 10^ 5 Hz one derives B ^ 6 x 10’ 3 t 3 ^ gauss. In a 

synchrotron model the two other Interesting parameters are Its size, R and the 
particle density, N Q such that N Q y" n dy « n(y) dy . We can take the standard 

3 

synchrotron ennssivity formula and derive the product R x n Q . For MK 501 we derive 
R 3 n Q -v 4.2 x 10 3 . In this type of model the loss times scale as v^; so, one 
might expect the X-ray to be more variable then the visible-light or UV. 

D. Other Possibil i ties 

1 . E mission from an Accretion Disk Around a Black Hole 

The spectral energy distribution in this case depends heavily upon a 
number of assumptions regarding the physical parameters involved. The calculations 
currently available (Eardley et al. 1978) are limited to a black-holes of mass 
100 M q and 10.0 M Q . The spectral turnover predicted by these authors for their 
standard disk and two-temperature disk models typically occurs in the far ultra- 
violet, while an optically thin model predicts no spectral turnover. Comparing the 
calculations for = 100 M Q and 10 M Q shows that the turnover point shifts 

toward longer wavelengths as the mass of the black-hole is increased, suggesting 
the possibility that for sufficiently large mass, the turnover may occur in ultra- 
violet or visible light. In order to account for the radio flux, however, another 
component (e.g., synchrotron emission) is needed in addition to the bremsstrahlung 
component. Further theoretical calculations are needed to extend the model to 
include masses for which the turnover occur in optical frequencies. 

2. E mission from Acceleration Limited Radiation 

Recently Cavaliere and Morrison (1980) have calculated the emission from 
a source where continuous reacceleration is taking place. In this model, visible 
light and UV emission can be characterized by F = A " a where a = (m-l)/(m-a), in 
which the magnetic field B has a radial dependence B = B Q (r/R)" ,n and P(r) = 

P Q (r/R)~ a is the radial dependence of the energy gain of a particle. There are 
3 possible higher energy components of slopes v = l/(. 9 m-a), or l/(?-a), or 1/a 
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depending on the dominant region In the model. If we let a be set by the radio 
emission and 6 by the X-ray we derive m ■ 1 , a ■ 1 , for all 3 cases. Thus the 
X-ray, UV and radio spectra uniquely define the radial dependences In this model. 

V. CONCLUDING REMARKS 

We have presented the first "quasl-lnstantaneous" spectral energy profile of 
the BL Lac source MK 501 and have compared It with predictions of two principal 
candidate emission mechanisms. 

While the synchrotron self-Compton model Is entirely consistent with our 
data and, as such, Is a favored model, an unfortunate result of synchrotron self- 
Compton theory Is that tight constraints on parameters cannot be achieved, as 
Is evident from Inspection of the relationships given in section IV (A). 

A synchrotron origin for the optical -X-ray spectrum also might be 
reasonable. 

Compton emission from thermal electrons can fit our data In the X-ray, 
ultraviolet and visible light regions, but we do have to assume an additional 
mechanism to account fur the radio flux. The required size light year) can 
be accommodated without necessarily assuming an exotic explanation such as a 
black hole. However, the source of soft photons needed for the Compton emission 
must also be explained. 

We also do not rule out other mechanisms, although direct comparisons 
are presently not practical. In particular, an accretion disk surrounding a 
massive black hole giving rise to bremsstrahl ung emission is a possibility to 
be further explored with theoretical calculations. 

Since we are dealing with a variable source in all frequencies, with the 
possible exception of the radio, we may unravel the energy mechanism by again 
observing the source simultaneously in all frequencies when it is at a different 
flux level. Any change In the spectral energy profile will provide additional 
eludes for further constraining the model. 

Since the UV and X-ray spectra imply continuity through the extreme UV to 
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soft X-ray, formal Integration of this spectrum shows that most of the energy 

emitted by MK 501 appears In the extreme UV. This suggests that the previous 

estimates of the total luminosity has been underestimated by a factor of about 

three to four resulting in an Integrated luminosity over our entire observed 

H “3 

bandwidth of ^ 1 x 10^ (g^) erg s" 1 . 

The HEA0-A2 experiment Is a collaborative effort led by E. Boldt of 
GSFC and G. Garmlre of CIT with collaborators at GSFC, CIT, JPL and UCB. Thanks 
are due Dr. A. Boggess and the IUE Observatory staff at GSFC for their competent 
assistance In obtaining the data. Y.K., R.F.M. , R.L.H. and K.R.H.H. are guest 
observers on the IUE. R.L.H. and K.R.H.H. gratefully acknowledge the assistance 
of W. Klnzel In the reduction of the IUE data and funding by NASA Grant NSG 5238 
and the WKU Faculty Research Grant Committee. G.N. and K.M. acknowledge support 
by NASA. We thank A.D. Code for helpful discussions of the interstellar reddening 
and E.A. Boldt for valuable comments on the X-ray results. P.A.F. wishes to 
thank C.R. Purton, P.C. Gregory, and W.M. Turner for assistance with the obser- 
vations at A.R.O. 
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TABLE 1. INFORMATION ON OBSERVATIONS 


Spectral Range 

log v 

Other Units 

Dates (1978) 

X-ray 

19.08-17.68 

50-2 keV 

Aug. 21-31 


19.08-17.68 

50-2 keV 

Sep. 8, 9 


uv 

15.36-15.20 

0 

1275-1910 A 

Aug. 

25 


15.20-14.97 

0 

1900-3140 A 

Aug. 

26 

Visible Light 

14.96-14.42 

0 

3290-11410 A 

Aug. 

29, 


IR 

1-4.13 

2.2 mu 

Aug. 

31 

Radio 

10.03 

10650 MHz 

Aug. 

30 - Sep. 3 


10.02 

10275 MHz 

Aug. 

15-24 


9.68 

4750 MHz 

Aug. 

15-31 


instrument 
HEAO-1 A2 
HEAO-1 A2 


IUE 

IUE 


5-meter 


2.5-meter 


46-meter 

46-meter 


300- ft 


FIGURE CAP T IONS 


Figure 1 - The X-ray spectrum of MK 501 In 1978 Aug-Sept, The upper left 

panel shows the best fit power law spectrum folded through the detector 

response compared to the pulse height data The lower left panel gives 

the incident spectrum unfolded using the best fit model. The lower right 

panel shows the 90% confidence contour for the X-ray power law Index 

?1 ? 

and the X-ray absorption In units of 10 at/cm „ 

Figure 2 - Ultraviolet spectral energy distribution of MK 501 from IUE obser- 
vations on 1978 August 25-26 Open circles represent data without correction 
for reddening. Solid circles indicate data corrected for reddening 
assuming E(B-V) = 0 083 mag. The points represent averages for frequency 
bins containing 10 points each, with error bars indicating + 2o Internal 
error for the bin means. 

Figure 3 - Absolute spectral energy distribution of MK 501 in visible light 
on 1978 August 30 through an aperture of diameter 10" on the 5-meter 
telescope. 

Figure 4 - 46-meter radio observations of MK 501 at 10650 (a) and 10275 (0) 

MHz obtained on 1978 August 15 through September 3. Error bars on individual 
observations are + 2o. The dashed line indicates the mean flux density, 

1 o 14 Jy. 

Figure 5 - 92 -meter radio observations of MK 501 at 4750 MHz obtained on 
1978 August 14 through September 1. Error bars are + 2 o„ 

Figure 6 - Combined absolute spectral energy distribution of MK 501 from 
X-ray to radio wavelengths based on the observations obtained in 
1978 August - September, 
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